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Thermodynamics of a heavy quark-antiquark pair in SU(3)-QCD is studied below the deconfinement critical 
temperature, T c . In the quenched case, a model of the string passing through heavy valence gluons yields a correct 
estimate of T c and a behavior of the string tension near T c . For two light flavors, entropy and internal energy 
can be obtained from the partition function of heavy-light mesons and baryons. They are in a good qualitative 
agreement with the lattice results. 



1. INTRODUCTION 

In the last few years, various indications ap- 
peared that the quark-gluon plasma is not a 
weakly coupled system even up to temperatures 
of the order of a few times the deconfinement one, 
T c . For instance, RHIC data point to a nearly 
perfect-fluid type behavior of the plasma at such 
temperatures, which is characterized by the shear 
viscosity to the entropy-density ratio of the order 
of 0.1, whereas the values for this ratio predicted 
by perturbative QCD are larger than 1. 

Besides experimental data, a lot of information 
on the nonperturbative properties of the plasma 
comes from the lattice. Among the recent results 
obtained there and calling for a theoretical ex- 
planation, are the anomalously large (from the 
perturbation theory standpoint) maxima of the 
entropy and internal energy of the static quark- 
antiquark pair in unquenched QCD around T c pQ, 
0. This work is aimed at a description of these 
data below T c . The strategy of our analysis is 
the following. In the next section, we will deter- 
mine T c and the effective string tension <r(T) in 
the quenched SU(A C ) QCD within the so-called 
gluon-chain model. In section 3, with the use of 
cr(T) adjusted to the unquenched (N c = 3, Nf = 
2)-case, we will calculate, within the relativistic 
quark model, the partition function of heavy- light 
mesons and baryons. Entropy and internal en- 
ergy stemming from this partition function will 
further be compared to the corresponding lattice 



data, and in the Summary section conclusions will 
be drawn. 

2. STATIC QQ-PAIR WITH GLUONS 

It is generally accepted that the QQ-string, 
which sweeps out the flat surface of the cor- 
responding Wilson loop, is produced by soft 
(stochastic background) gluonic fields. On top of 
those, fluctuations of the gauge field exist, which 
lead to string vibrations. These fluctuations are 
called valence gluons, and the QQ-stfrng pass- 
ing through one such gluon is called hybrid. At 
asymptotically large QQ-separations of interest, 
one has a string passing through numerous va- 
lence gluons, which we will call gluon chain. Fur- 
ther, when one starts heating the system, at low 
temperatures the free energy of one string bit be- 
tween two nearest gluons in the chain, being a 
constant, exceeds thermal mass of a gluon, which 
grows with T linearly from zero on. The motion 
of gluons at low temperatures is therefore guided 
by the string - they move collectively and do not 
affect string's dynamics. However, at a certain 
temperature To smaller than T c , the gluon's ther- 
mal mass becomes larger than the free energy of 
one string bit. From this temperature on, the 
system looks totally different, since gluons from 
the string's standpoint are now (nearly) static. 
Therefore, at To < T < T c , a gluon chain is noth- 
ing but a sequence of static nodes with adjoint 
charges, connected by independently fluctuating 
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string bits. At the moment of formation of such a 
chain, its end-point originating from the heavy Q 
performs a random walk towards Q over the lat- 
tice of static nodes. The entropy of such a random 
walk turns out to be large, namely proportional to 
its length, and eventually leads to the deconfine- 
ment phase transition in this model. The reason 
is that color may alter from one node to another, 
i.e. every string bit may transport each of the N c 
colors. The total number of states of the gluon 
chain is therefore N^ a ', where L is its length and 
a is the length of one bit. 

Taking into account that the full free energy 
is the sum of the usual linear potential and the 
free energy of the random walk, we have for the 
effective string tension: 
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is the partition function of the random walk. 
Here, j3 = 1/T, a = (440 MeV) 2 is the zero- 
temperature string tension, s — aL is the 
Schwinger proper time, and n is the number of 
a Matsubara mode. This formula takes into ac- 
count the confining potential between the end- 
point of the walk and its starting point (chosen 
at the origin), as well as the entropy factor dis- 
cussed above. At asymptotically large R's of in- 
terest, only the contribution of the (n = 0)-term 
is essential, and we obtain 
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An estimate for T c stems from the condition that 
the argument of the first square root vanishes: 
era 

w c >i ln7V c ' 
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Equating T c to the modern N c — 3 lattice 
value 1 , 270 MeV, we obtain an effective length 
of one string bit a ~ 0.31 fm. This is larger 
than the minimal possible value of this quantity, 
a = 0.22 fm - the so-called vacuum correlation 
length H3, which defines the onset of a string- bit 
formation. The temperature To, below which the 
lattice of valence gluons does not exist, can be de- 
fined from the condition cr(T c ) — 0, which yields 
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130 MeV. 



An important finding of this model is the behav- 
ior 



tr(T) - ^Tc-T at T c . 

It is the same as the one which follows from the 
Nambu-Goto model for the two-point correlation 
function of Polyakov loops 0] . 

Let us finally consider the limiting case when 
string bits cannot alter color, i.e. one should for- 
mally set in Eq. N c = 1, that yields 
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Determining in this case To directly from the 
equality of the gluon's thermal mass in QCD to 
the free energy of one string bit, we have 
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Setting for an estimate a = 0.22 fm, g = 2.5, we 
obtain still reasonable values T c = 290 MeV, To = 
85 MeV. However, the fundamental difference of 
this limiting case from the realistic one, N c > 1, is 

that herea(T) -» \ yj a(1+ a ^ ] (T C ~T) ~ (T c -T) 
at T — > T c . Such type of the critical behav- 
ior characterizes the 2d Ising model, which by 
no means can be realized in the 4d quenched 
SU(Ay QCD. The same linear fall-off of a(T) 
with (T c — T) one finds also in the Hagedorn 
phase transition and in the deconfinemcnt sce- 
nario based on the condensation of long closed 
strings [SJ. It is a mere consequence of the for- 
mula cr(T) = a - 2j£ at S oc R. 
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3. STATIC QQ-PAIR 
QUARKS 



WITH LIGHT 



Let us now consider the unquenched case, 
where, at a certain distance, the QQ-string breaks 
due to the production of a light qq-pair. The sub- 
sequent hadronization process leads to the for- 
mation of heavy-light mesons {Qq) and heavy- 
light- light baryons (Qqq), as well as their an- 
tiparticles. We will consider the (Nf — 2)-case, 
with light u- and (i-quarks, and use the value 
T c = 200 MeV 0. Equation © then yields an 
effective value of a = 0.23 fm, which should now 
be used in Eq. to determine cr(T). Then, for 
a heavy-light meson, the Hamiltonian of the rel- 
ativistic quark model reads 



H Qq 



Q 



V(r). 
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Here, V(r) = a(T)r - (2 - 5)^/^{T), m Q is the 
heavy-antiquark mass, whereas m q is the con- 
stituent mass of a light quark, m q ~ 300 MeV. 
The subtraction of 2^J a(T) in V(r) is known to 
be important to reach an agreement between the 
predictions of the relativistic quark model to the 
phenomenology of meson spectroscopy [H| , [3 ■ An 
additional correction, 5^J <x(T), is needed to obey 
the normalization condition F — » 2(mjjo —m c ) at 
T — > 0. Here, F is the resulting free energy of two 
noninteracting mesons, m D o — 1.864 GeV is the 
£>°-meson mass, m c = 1.48 GeV is the c-quark 
mass. Further, the square root in the relativistic 
kinetic energy can be disentangled by using the 
integration over an auxiliary parameter: 
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One then arrives at a 3d Schrodinger equation of 
the form 

(Kd 2 + M |x|)i/»(x) = £ty(x), 

where K and M are positive constants of di- 
mension [length] and [mass] 2 , respectively. Its 
eigenenergies have the form \7\ 



E nr ,i = a nr ,i {KM 1 ) 



where a nr /s are positive numbers with a 00 ~ 
2.34. Further strategy for the calculation of the 
partition function of a meson \J\ consists of ac- 
counting for its ground state exactly and model- 
ing its higher eigenenergies by those of a 3d har- 
monic oscillator with the frequency 



lo = {KM 2 ) 1/3 



4^ 2 



1/3 



Then, because of the degeneracy factor (n/2 
l)(n + 1) for the oscillator, we have 

Z Qq = -2= exp (-(3m Q + Ip^T)) x 
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The free energy of two noninteracting mesons 
reads F — ~ Tln(4Zg ), where the factor 4 is a 
product of two flavors and two spin states of the 
mesons (which are antiparticles of each other). 
Doing the sum over n analytically, we arrive at 
the following result: 
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The remaining /i-integration has been done nu- 
merically. The resulting value of the parame- 
ter 8 is 0.344. Entropy and internal energy can 
further be calculated by the standard formulae 
S = -dF/dT, U = F + TS. Note that, when cal- 
culating the entropy, the derivative d/dT should 
not act on the Hamiltonian in the partition func- 
tion. 

Similarly, one can treat (Qgq)-baryons, whose 
Hamiltonian reads 
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Here, we have approximated the position of the 
baryon string-junction point by the position of 
Q, that is legitimate due to the heaviness of Q. 
The value of the parameter S, at which the free 
energy of a baryon goes to m(A+) — m c at T — > 0, 
is 0.393. Here, m(A+) = 2.286 GeV is the mass of 
the -baryon. In the figures, we plot the entropy 
and the internal energy of mesons together with 
baryons, at the averaged value S — 0.37. In the 
same figures, we plot the corresponding lattice 
data for the QQ-pair FJ]. 

4. SUMMARY 

In this work, we have analytically addressed 
thermodynamics of the static QQ-pair below 
the deconfinement phase transition, in quenched 
and unquenched QCD. In the quenched SU(_/V C )- 



case, within the so-called gluon-chain model, we 
have derived an effective temperature-dependent 
string tension. With this string tension adjusted 
to the unquenched (N c — 3, Nf = 2)-case, we 
have calculated the entropies and internal ener- 
gies of heavy- light mesons and baryons. In gen- 
eral, the results obtained are in a good agree- 
ment with the corresponding recent lattice data 
for the same quantities. Our analysis shows also 
that the entropy and the internal energy of heavy- 
light mesons alone are not enough to reproduce 
correctly the lattice data. 
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